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In this paper, we focus on the modeling and control of a wind power system 
based on a Permanent Magnet Synchronous Generator (PMSG). We 
proposed a technique of control strategies to have the maximum power from 
wind turbine (WT). This study deals with the problem of Maximum Power 
Point Tracking (MPPT) based on Takagi Sugeno fuzzy model. The stability 
analysis is achieved. The gains of the designed controller are calculated by 
solving Linear Matrix Inequality (LMI). Finally, simulation results are 
provided to demonstrate the validity and the effectiveness of the proposed 
method. 
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1. INTRODUCTION 

A serious damage of the environment is caused by global warning and fossil fuel pollution and in 
consequence it became a problem for humanity. These reasons led researchers to present alternative sources 
of supply such as photovoltaic and wind energy systems. These solutions are considered as the most 
sustainable sources for many uses for example irrigation and electrification in isolated areas [1]. Also, 
Variable Speed Wind Energy Conversion Systems are the dominant technologies in the present wind power 
industry because they have multiples advantages, over the fixed velocity systems, as the ability to obtain 
Maximum Power Point Tracking (MPPT) control methodology in order to extract maximum power at 
different wind, higher overall efficiency, power quality and it can be controlled to reduce aerodynamic noise 
and mechanical stress on VS-WECS by absorbing the wind-power fluctuations. Moreover, the Wind 
Turbines Generators based on Permanent Magnet Synchronous Generators (PMSG) are considered popular 
for variable-speed generation system and the use of the PMSG in large WTG is increasing. It is connected 
directly to the turbine without gearbox and so it can operate at low speeds [2]. In addition, it can decrease 
again weight, losses, costs, demands maintenance requirements. 

In order to extract the optimal accessible power given by the WECS, several studies have presented 
multiple control schemes [3]. The simplest techniques for searching the maximum power point are based on 
PI controller but this method is still classical and lack of performance, the perturb & observe method is the 
simplest algorithm for searching the maximum power point but the duty of the power converter perturbs 
rapidly and causes energy loss near the maximum power point [4]-[5], The neural network method estimates 
the wind speed from the measured turbine power and generator speed and thus decides the maximum power 
generator speed command or torque reference for the operational point tracker [6]. In addition, fuzzy logic 
methods have been widely adopted and they have presented better effectiveness [7]. Nevertheless, the main 
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issues of most of these methods are lack of stability and strict theoretical analysis so that the maximum power 
point varies over a wide range. 

In literature, the T-S fuzzy model-based control has been fully mentioned for example [8]-[9]. Thus, 
the control based on (T-S) model has been more popular as one of the most successful techniques for systems 
and control applications due to its reliability and effectiveness. The principal advantage is that the controller 
is systematically realized by using parallel distributed compensation (PDC) and a linear matrix inequality 
(LMI) technique [10]. The fuzzy model of Takagi-Sugeno allows modelling the nonlinear system by 
representing local dynamics by linear models. As a result, the global system model is achieved by a 
combination of the different linear models. Then, a linear feedback control has to be constructed for every 
local linear model. Accordingly, the consequent overall nonlinear controller is once more a fuzzy 
combination of each distinct linear controller [11]. 

The rest of this paper is arranged as following: in the second section, a dynamic modelling of the 
WECS is presented. The third section starts with presenting the control strategy. Then, fuzzy modelling of 
the system is presented. Finally, the control design is explained and stability conditions are given. Section 4 
is reserved to present numerical simulations that illustrate the effectiveness of the proposed control topology. 
We finish by a conclusion. 


2. MODELLING OF THE WIND TURBINE GENERATOR 
3.1. Wind turbine model 

A schematic overview of the WECS is shown in Figure 1. The system supplies a resistive load and 
consists of a wind turbine rotor, PMSG, rectifier, and a boost converter. Wind turbine converts the wind 
energy into mechanical energy, which then runs a generator to create electrical energy. 


Turbine 



Figure 1. Structural diagram of WECS 


The amount of power captured by the wind turbine is given as [12]: 

P m = X -pnR^C p {A) ( 1 ) 

Where, P is the air density ( kg/m3 ), R is the blade radius (m), V is the wind velocity (in m/s), and C p (A) 

is the power coefficient usually given as a function of the tip-speed ratio A .The tip-speed ratio k is defined 
as 


A = 


Ro \, 

V 


The wind turbine mechanical torque output Tm given as: 


l 5 c„U) 2 

1 = — P7TK ——m 

jti 2 ' ^£,3 ^ 


( 2 ) 


( 3 ) 


Many different versions of fitted equations for C p have been used in previous studies. This paper defined 
C p based on the following [13]: 
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Cp = -0.212/t 3 +0.0856/1 2 +0.2539/1 (4) 

The maximum of C p , that is C Pma =0.15 , is reached for X opt =0.78. Hence, there is one particular 
X = X opt , optimal power coefficient, and an optimal turbine rotating velocity CO m opt , under which C p takes 
a maximum value C p max . X opt is the optimal value of the tip speed ratio, which is dependent on the 

characteristics of the turbine system. Also, the maximum power can be captured from the wind. Accordingly, 
the turbine system works in maximum power point tracking (MPPT) [14], as depicted in Figure 2. 

The Figure 3 presents C p as a function of X . According to the figure, there is only one optimal 

point, denoted byX t , where C p is maximum. Continuous operation of the wind turbine at this point 

guarantees that it will obtain the maximum available power from the wind at any speed, as shown in Figure 
4. 



Generator Speed (rad/sec) 



Figure 2. Characteristics of turbine power as a 
function of the rotor speed 


Figure 3. The characteristic of the power coefficient 
as a function of tip speed ratio 


3.2. PMSG Model 

In order to simplify the dynamic model of the electrical part of PMSG, applying a change of 
landmarks which rotates at the same speed as the rotor (Two bobbins dq rotate with the rotor and generating 
the same effect as the three fixed coils) [14]. Then the electrical model of the PMSG is defined in the rotating 
reference fame d-q as follows: 

dI " -co L 1 (5) 

p n an v ' 


V sd - R Jsd + L d 


sd 


dt 


V sq =R s I sq +L q 


dl 


dt 


- + CO e( L d I sd + y / m) 


( 6 ) 


Where V sq , V sd are the direct and quadrature components of the PMSG voltages, R s , L d and L , 
respectively, are the resistance, the direct and the quadrature inductance of the PMSG winding, y/ m (wb) 
represents the magnet flux, CO e (rad/s) is the electrical rotational speed of PMSG, I sd , I sq (A) are the direct 
and quadrature components of the PMSG currents, respectively. 

The dynamic equation of the wind turbine is described by: 


m, 

dt 


=—t +^n t — 1 -t, 

J T J T J T 

(7) 
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The electromagnetic torque of a p-pole machine is obtained as [16]: 

T =-n (w I +(L,-L)I,I ) (8) 

Where T em (N m) is the electromagnetic torque, n p is the number of pole pairs, D T is the damping 
coefficient and J T is the moment of inertia. 

3. CONTROL STRATEGY 

3.1 The control strategy of the MPPT 


Turbine 



Figure 4. Fuzzy tracking control diagram 


The Figure 4 present the control scheme treated in the current work in the aim to maintain generated 
power at its maximum value. MPPT is a process that searches the operating point which allows extracting the 
maximum available produced power at any wind speed level. The control algorithm requires to measure 
reference signals to be tracked for the rotor speed, the direct and quadrature stator current. Then, the MPPT 

control is realized by acting on the rectifier side by generating control signals ( U dr and U qr ) to track the 
peak power point. The rotor speed is measured according to the following equation, 


Q 


r 


4A 

R 


(9) 


In this study, the d-axis reference current is chosen as: i sdr = 0 we can calculate the reference of the current 
in quadrature given, once the rotor speed reference isdetermined by: 


C 


w 


e -^L = Cte.^ 


Where 


k,=(p f xp 


And 


Cte = 


pxR, C 


/?max 


2/t. 


opt 


( 10 ) 


3.2 T-S Fuzzy model of wind energy conversion system 

In the T-S fuzzy models, the nonlinear system of a SPMSM is represented by several linear sub¬ 
systems according to the model rules. Each sub-model leads to the overall behavior of the nonlinear system 
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using a weighting function [2]. The dynamic model of the PMSG-WT can be described by the following 
nonlinear state space form [13]: 


x(t) = A(x(t))x(t) + Bu(t) + ET m (t) 


( 11 ) 


Where: x — ( Q / , i ) 

\ m sd sq J 


( 12 ) 


x(t ) are the state variables, u(t) are the control inputs, A, B i et E { are the state matrices of the sub-system. 
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Where V d and V represent the average control signals of the inverter in d-q frame. Moreover, the 

nonlinear model (12) can be presented by a T-S fuzzy model with r = 2 1 fuzzy If-Then rules as 
follows: 


Model rules i: 

If z(t) is F/ then x(t) = Ax l (t) + B^uit) + E^T m {t) (13) 

If zit) is F\ then x(t) = Ax 2 (t) + B 2 u(t) + E 2 T m (t) (14) 


Wherez(f) = CO m (t)is the premise variable and F 1 is the membership functions, Let define 
D = m&xiojjt)) and d = min(co m (t)) , so CO m (t) = [d,D ]. 


1 D-d 
F^l-F/ 

Hence, the local subsystem matrices are given by: 
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(15) 

(16) 
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After defuzzification, the fuzzy system of a PMSM can be expressed as: 


kt ) = 2>,.(z(0)(4x(0+ B M 0 + ETjty) 

i =1 


(17) 


Where 


h,(z(t)) = 


F]pm 

±F;(z(t)) 

7=1 


Vf > 0, hj ( zit )) > 0 and ^ h t (z(t)) = 1 

i =1 


(18) 


3.3 The PDC control 

The goal is to drive the actual variables X of PMSM-WT system to track the optimal trajectory X opt , such 
that: 

(0 = x(t) -X opt (0 —>0as ? —>00 (19) 


Where x(t ) = CO m (t) , x t (£) is the desired value of , and X, (7) be defined as the tracking error and its 

time derivative is given by: 

• • • 

( t) = x(t)~ Xop it) (20) 

Hence, the new PDC fuzzy controller is designed such as [13]: 

2 

T = ~Y J h i (z(t))K i x e (t) (21) 

1=1 

and r = [r q ,T d J (22) 

Where T q , T d are new controller to be designed via LMIs approach, K { denotes the control gain 

corresponding to each linear sub-model and .£,(/) is the augmented state vector. Therefore, by replacing T 
by its expression, the closed loop model is described by [15]: 


x e (t) = Hh i (z(t))h j (z(t))(A i -B^^x^t) 

i =1 7=1 

With: G^iA-B'Kq) 

• J _ r 

The can be written as: X e {t) — 

i =1 7=1 


(23) 


(24) 
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4. RESULTS AND DISCUSSION 

In this section, a numerical simulation was carried out by using MATLAB/Simulink, the system 
described in Figure 4. First, we have simulated the system without controllers. The wind speed is given as a 
sum of several harmonics [16]: 

n 

Kent =V 0+ 5X sin (fflL .t) , (Figure 5). (25) 



Figure 5. Wind speed variation 


The proposed fuzzy controller results in the MPPT responses shown in Figure 6, where the 
trajectory of the speed of the PMSM is tracking its reference signal and the Power coefficient is quickly 
achieved. Next, to carry out comparisons with the traditional control method, we apply proportional-integral 
(PI) control to track the Rotor speed and the Power coefficient Cp 

After applying the proposed fuzzy controller, the responses of the Rotor speed and the Power 
coefficient are shown in Figure 7, Figure 8.the control result is also compared with the PI based MPPT, 
where the comparison of the tracking error demonstrate that the system with the synthesized T-S fuzzy 
controller has a good behavior. Indeed, the rotor speed and the power coefficient track well the reference 
trajectory with good reliability over the whole speed range. 



Figure 6a. Power coefficient Cp Figure 6b. Rotor speed tracking 

Figure 6. Nonlinear T-S Fuzzy control method 
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Figure 7a. Power coefficient Cp of PI and TS methods 


Figure 7b. Zoom of coefficient Cp 


Figure 7. Response of Power coefficient Cp 



Figure 8a. Rotor speed tracking of PI and TS methods 



Figure 8b. Zoom of rotor speed tracking 


Figure 8. Response of rotor speed tracking 


5. CONCLUSION 

This work has been devoted to modeling, simulation and analysis of a wind conversion system with 
the main objective to ensure maximum peak power tracking (MPPT). We have proposed a fuzzy tracking 
control for a PMSG-WT, the proposed controller has been compared with the PI controller. Simulation 
results show that the proposed wind system gives best performance and were found the faster than the PI 
controller. 
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